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PSEUDOCHAOS 

G.M. Zaslavsky 
M. Edelman 



ABSTRACT A family of the billiard-type systems with zero Lyapunov 
exponent is considered as an example of dynamics which is between the 
regular one and chaotic mixing. This type of dynamics is called "pseu- 
dochaos". We demonstrate how the fractional kinetic equation can be in- 
troduced for the pseudochaos and how the main critical exponents of the 
fractional kinetics can be evaluated from the dynamics. Problems related 
to pseudochaos are discussed: Poincare recurrences, continued fractions, 
log-periodicity, rhombic billiards, and others. Pseudochaotic dynamics and 
fractional kinetics can be applied to streamlines or magnetic field lines be- 
havior. 



1 Introduction 



There are many examples that show transition to turbulence through the 



spatio-temporal chaos. The case of 2-dimcnsional Kolmogorov flow (Piatt 



3irovich, and Fitzmauricc [1991 1) displays a complicated alternation of dif- 
ferent laminar and chaotic regimes following changes of a control param- 
eter, which precede the dynamics known as the turbulent one. While the 
common opinion focuses on an exceptional role of coherent structures pre- 
sented even in turbulent flows, the structures are still ellusive and we still 
do not have a universal method for their depiction. Interesting ideas of 
t he diagnostic o f coherent structures can be found in recent publication s 
( |5irovich [198j , |Weiss, Provcnzale, and McWilliams [1998| , |Hallcr [200(1 ). 
An indirect analysis of flows can be developed from the trajectories of pas- 
sive particles (tracers). This aproach is known also as the Lagrangian one 
versus the Eulerian one. 

While the study of tracer dynamics is not sufficient to provide full in- 
formation about the corresponding flow, it can reveal an important in- 
formation about the flow transitions from one regime to another. An im- 
portant example of the tracer dynamics appears for the so-called ABC 
(Arnold-Beltrami-Childress) flow ( Arnold [1978 1) for which the tracer dy 



namics coincides with streamline trajectories, is Hamiltonian (Zaslavsky, 



3agdeev, Usikov, and Chcrnikov [1991]), and consists of regular and chaotic 



domains in phase space for all nonzero values of the flow parameters (Za 



lavsky, Sagdeev, Usikov, and Chcrnikov [1991], Hcnon [1966]). An area 
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of chaotic streamlines depends on the symmetry of flow, flow parameters 
spectra, etc. ( Agullo, Vcrga, and Zaslavsky [1997|, Benkadda, Gabbai, and 



Zaslavsky [1997|7Benkadda and Beyer [2002[ |, Zaslavsky, Sagdccv, Usikov 
and Chcrnikov [f99fp ~ 



Similar to the tracer dynamics problem is the behavior of magnetic field 



lines for a given vector-field B(r) ( Roscnbluth, Sagdeev, Taylor, and Za- 
slavsky [1966 1 , jZaslavsky, Sagdccv, Usikov, and Chcrnikov [1991| ) known 



to be chaotic long before the studies of chaos of Lagrangian particles. This 
analogy has been used to study a novel type of problem: randomness with- 



out sensitiv ity to the perturbation of initial conditions ( Zaslavsky and Edcl 



man [2001 ] ) . Field- lines of v or B can wind invariant surfaces of a topolog- 



ical genus more than one. In this case equations for the field lines are not 
integrable even if the Lyapunov exponent ctl = ( Kozlov [1996 1 ) . We call 
trajectories to be pseudochaotic if <jl = but the equations that define the 
trajectories are not integrable. 



It w as shown for some cases of the pseudochaos ( Zaslavsky and Edclman 



[2001 1) that the ensemble of trajectories can be considered in a statistical 



manner and can be described by a kind of kinetic equation with fractional 
derivatives. Kinetics of trajectories reflects self-similarity of dynamics and 
supcrdiffusivc equations for the moments of a distribution function along 
coordinates, i.e. supcrdiffusivc transport. This article continues the work 



Zaslavsky and Edclman [2001] replacing the problem of the field- lines by a 
billiard type problem. There is no one-to-one correspondence between the 
field-line trajectories problem and the billiard one. Nevertheless, the study 
of particle dynamics in the billiard-type systems can provide a realistic 
insight on different possibilities of the dynamics as it has happened with 
the Sinai billiard. 

We consider different billiard- type models with zero Lyapunov exponent. 
These models can be applied straightforwardly to such physi cal problems as 



ray dynamics in a media with complex (fractal) boundary ( Sapoval, Gob 



ron, and Margolina [1991 ], Hebcrt, Sapoval, and Russ [1999| ) , light (Wilkin- 



son, Fromhold, Taylor, and Micolich [2001a , b|) and sound ( Zaslavsky and 



Abdullacv [1997 ] ) propagation in nonuniform med ia, magnetic field lines in 
toroidal plasmas ( ^aslavsky and Edclman [2001 1), tracer dynamics in re- 
connected vortices ( Zaslavsky and Edclman [2001 1), et c. Similar problems 



prove d to be interestin g in problems of quantum chaos (Richcns and Berry 



[198l[, |Wicrsig [2000|, |Artuso, Guarncri, and Rcbuzzini [2000| , |Artuso 



Casati, and Guarncri [1997 1 ) and general p roblems of anomalous trans 
port ( Artuso, Guarncri, and Rcbuzzini [2000 1 , Artuso, Casati, and Guarncri 



1997| , |Zwanzig [1983| , |Zaslavsky and Edclman [2001 



portant works on the polygonal billiards ( Galpcrin and Zcmlyakov [199C | 



There are also im- 



Katok [1980(| , [Katok and Hassclblatt [1995| , |Gutkin [1996|) and related 



to it. The interval exchange transformation (IET) (Katok [1980 1 , Zorich 



[1997 1 ) which considers complexity, mixing, and transport. 
In this paper we briefly review our previous results on the kinetics and 
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transport properties of the trajectories in a family of polygonal billiards 
and the corresponding Lorentz-type gases ( Zaslavsky and Edchnan [2001[ |), 
and present new results on the transport exponents and more complicated 
billiard-type models. It is worthwhile to mention here, that while there 
is a way of evaluating the transport exponents for IET ( Zorich [1997| ), 
these results cannot be applied immediately to the transport exponent of 
trajectories, and a special analysis of trajectories and their ensemble is 
necessary. 

There are three separate problems related to the pseudochaotic kinet- 
ics, that will be briefly discussed here: transport in the rhombic billiard, 
log-periodic properties of continued fractions, and the origin of specific ex- 
ponents for Poincare recurrences. 



FILAMENTED AND MERGED 
SURFACES, BILLIARDS, AND WEAK 
MIXING 



We will speak about the particle dynamics having in mind a tracer or a field 
line as a particle trajectory. Examples of the filamentcd and merged surfaces 
are given in FIGURES 2.1(a) and (b) correspondingly. The case (a) appears 



in the toroidal palsma confinement (Zaslavsky and Edchnan [2001 1, Moro 



zov and Solovicv [1966]) while the case (b) can appear in both plasma and 
fluids. Under some conditions, the problem of the geodesies along surfaces 



can be reduced to the problem of trajectories in cllastic billiards (Galperin 
and Zcmlyakov [1990 1, |Katok and Hassclblatt [1995 1 ) . Some examples of 
billiards are given in FIGURE 2.2 Billiards must have rational angles to be 
reducible to the equivalent dynamics along surfaces. Each of the billiards in 
FIGURE 2.2 can be periodically continued in two or one directions (FIG- 
URES 2.3 and 2.4) creating a "generalized Lorentz gas" (GLG). As an 
example, let us point an evident connection between geodesies in FIGURE 
2.1(a), and trajectories in FIGURES 2.2(a) and 2.3(a). 

Billiard trajectories can be studied by introducing a map of some interval 
into itself. For the case in FIGURE 2.2(a) one can consider the map of the 
slit into itself, which belongs to the type of the interval exchange trans- 
formation: (IET) ( |Galpcrin and Zcmlyakov [1990H , |Zorich [1997| ). There 
are some common features between all four types of problems mentioned 
above: geodesies along complex surfaces, billiards, GLG, and IET. All of 
them correspond to the dynamics with zero Lyapunov exponent, i.e. two 
initially close trajectories diverge not faster than a power of time. All of 
them describe the so-called pseudo-integrable (but not intcgrablc) situation 
( [Katok [1980| , |Katok and Hassclblatt [1995j] ) . Al l of them correspond to 



the dynamics with the weak-mixing (Katok [1980 1 , Katok and Hassclblatt 
[1995 1 ) (more accurate formulation can be found in |Gutkin [1996| ]) 



4 
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Let us recall ( Cornfcld, Fomin, and Sinai [1982| ] ) that for almost arbitrary 



square integrable functions G\{x) and G 2 {x) of the coordinate x in phase 
space, the weak mixing means that 

n— 1 

lim - V[(Gi(f k x)G 2 {x)) - (G 1 )(G 2 (x))} 2 = (2.1) 



while the mixing, or strong mixing, means 

lim [(G 1 (f"x)G 2 (x)) - (G 1 )(G 2 (x))} = (2.2) 

n — >oo 

where T is a time-shift operator and (...) means the averaging over a 
corresponding measure. Weak mixing can be accompanied by arbitrary 
large and long-lasting bursts, i.e. fluctuations with large deviations from 
zero of the correlation function 

R n = (f n x ■ x) - (x) 2 . (2.3) 

These fluctuations exist simultaneously with the property that time average 
of is zero. A sample of trajectory for the billiard in FIGURE 2.2(a) or 
GLG in FIGURE 2.3(a) is given in FIGURE 2.5 indicating arbitrary long- 
lasting almost periodic pieces of the trajectory. 

There is no one-to-one correspondence between four different types of the 
above-mentioned systems: geodesies on a compact surface, billiard, GLG, 
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FIGURE 2.2. Examples of billiards of different types. 



and IET. For example, properties of the IET obtained in Zorich [1997 1 
are not sufficient to describe kinetics of trajectories in the corresponding 
billiard since the equation for the mapping time is not involved in IET. 
This explains a major difficulty in studying the billiards and GLG. Let 
us mention that I ET dynamics and polygonal billiard dynamics is ergodic 
and weak mixing ( Katok [1980 1, Gutkin [1996 1 ) , but there is only a conjec- 
ture that the complexity of the billiard trajectories is at most polynomial 
( Gutkin [1996(1 , Zaslavsky and Edclman [2001 1). Due to this, we will call 
pseudochaotic the dynamics of trajectories in billiards and GLG, when their 
Lyapunov exponent is zero and their complexity is, at most, algebraic. As 
it will be shown, such billiards can reveal fairly good statistical properties. 

Let us provide a few examples for the billiard wit h a slit in FIGURE 
2.2(a) and the corresponding GLG in FIGURE 2.3(a) Zaslavsky and Edcl 
man [2001 1. Consider an element £ of the slit in FIGURE 2.2(a) and intro 



ducc a normalized distribution of the Poincare recurrences to £ as 



P Tec (t) = lim -P re 



(2.4) 



where P lec (£,t) is a probability density to return a trajectory to the ele- 
ment £, does not matter from which side of £. Due to the above-mentioned 
conjecture, we can expect that 



P rcc (i) - 1/f , (t oo) 
with a recurrence exponent 7 that satisfies the condition 

2 < 7 



due to the Kac lemma ( Cornfeld, Fomin, and Sinai [1982 ]) 

p 00 

r rcc = / dt t P(t) < 00 , 
Jq 



(2.5) 



(2.6) 



(2.7) 
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FIGURE 2.3. Periodically continued billiards form a "generalized Lorentz gas" 
(GLG). 




FIGURE 2.4. Ray propagation billiard model with a piecewise periodic nonuni- 
formity. 



i.e. the mean recurrence time is finite for the bounded area preserving 
dynamics. 

To study transport properties in billiards, consider a bundle of trajec- 



tories, launched with initially close angles, as an ensemble (Zaslavsky and 



Edclman [2001 1). For the same case in FIGURE 2. 3(a) we can consider mo 
ments (y 2m ) with an integer m. It was shown in Zaslavsky and Edclman 



[200l| that 



(y 2m ) ~ f^( m ) 
with 

fx(m) ~ m/i(l) , w 1.7 

and with the connection 

7« M (1) + 1 



7 w 2.7 



(2.8) 
(2.9) 
(2.10) 



Similar value of /i(l) was obtained in Artuso, Guarncri, and Rcbuzzin: 



[2000 1 . Some deviations of fj.(m) from the linear law (2.9) will not be dis- 
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FIGURE 2.5. Samples of a trajectory with different space (time) scales for the 
billiard in FIGURE 2.2(a). 



cussed here (see more in Zaslavsky and Edclman [2001 1). The values (2.10) 
arc almost independent on the length of the slit and the launching an- 
gle. This universality will be commented later as well as deviations from 
(2.8)-(2.10). The results (2.8)-(2.10) can be immediately transferred to the 
billiards in FIGURE 2.2(b) and to the GLG in FIGURES 3(b) and 3(c), 
since trajectories do not change initial tangent of the angles during the 
scattering. Formula (2.8) shows anomalous transport (/a(1) ^ 1) of the su- 
perdiffusion type and an approximate self-similarity since fJ,(m) ~ m/i(l). 
All these properties and deviations from them will be discussed in the 
forthcoming sections. 



3 MORE BILLIARDS 

In this seciton we would like to present more different types of billiards 
and their statistical properties. These results can be considered as a kind 
of "experimental material" . 

Let us start from the two equal-squares-lattices in FIGURES 3(e,f) with 
infinite and finite horizons, respectively. Some squares are rotated by an 
angle of irrational tangent with respect to the others. Due to that, almost 
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all trajectories rotate ergodically in the coordinate space in contrary to 
the billiards in FIGURES 3(a-d). Samples of trajectories are given in FIG- 
URES 3.1, 3.2. There are two important features common for both figures. 
First, the trajectories have arbitrary long flights. We call "a flight" any 
long part of the trajectory that, after small scale averaging, is almost reg- 
ular. Such pieces of trajectories correspond to an intermittent dynamics 
and they are responsible for the distribution functions with tails and for 
fractional kinetics (Zaslavsky [1994a ,0, Saichcv and Zaslavsky [1997 1 ) . The 



trajectories have similarity in small and in large time and space scales. The 
self-similarity is less evident in FIGURES 3.1, 3.2 than in FIGURE 2.5. 
More delicate description of the trajectories can be obtained from simula- 
tions. 

In FIGURE 3.2 we present two types of statistical properties of tra- 
jectories for the finite horizon billiard. The moments of the coordinate 
displacements 

{x 2m ) » t^ {m) , (y 2m ) « t^ {m) (3.1) 

with 

Hx{m) = H y {m) = [i(m) f=a m/i(l) (3.2) 
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FIGURE 3.2. A trajectory with two, small and large, time intervals (two top 
plates), and the moments {x 2m ) (regular lines), (y 2m ) (dash-point lines), and 
distribution of the recurrences P(t). The data for the two bottom plates are 
obtained after averaging over 8,192 trajectories and time 3.4 ■ 10 s for each. Four 
curves of the moments vs. time correspond to m = 1, 2, 3, 4 from the bottom. 



and 

H{1) w 1.5 . (3.3) 

Considering one cell of the lattice (similar to the billiards in FIGURE 2.1), 
one can obtain a distribution function of Poincare recurrences (see (2.5)). 
The simulation gives 

7 w 2.4 (3.4) 

in a good agreement with the relation (2.10). Similar to (3.3), values of 
/i(l) ps 1.5 were obtained for the billiard with the infinite horizon (FIGURE 
2.3(e)) with the recurrence exponent 7 ps 2.5. 



4 CONTINUED FRACTIONS AND 
SCALINGS 



It is convenient to consider first a theory for the billiard with a slit in FIG- 
URES 2.2(a) or 2.3(a) since the kinetics in the billiard is one-dimensional 
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along y and velocity v x along x is constant. This billiard attracts the at- 
tention fairly long ago (see Zwanzig [1983|] and fol lowing publicatio ns [Han- 



nay and McCraw [1990fl, |Richens and Berry [198~l|], |Wiersig [2000[ , [Artuso 



Guarneri, and Rcbuzzini [2000|, Artuso, Casati, and Guarncri [1997]). In 



our work Zaslavsky and Edelman [2001] we use the properties of the con- 
tinued fractions ( Khinchin [19641 7 Levy [1937 1) to obtain the recurrences 
distribution and kinetic properties of the trajectories. In this sect ion we 
will repeat and extend the results of Zaslavsky and Edelman [2001 1 . 

Let us recall the samples of the trajectory in FIGURE 2.5. Trajectories 
with an angle i? to x will be called rational/irrational if tani? is ratio- 
nal/irrational. We will consider only irrational ensembles, i.e. bunches of 
trajectories with irrational t&w&j and "dj £ (?? — Ai?/2,^ + Ai9/2) within 
an interval A$. 

Let 

tani9 = o + £ , (4.1) 

where ao is the integer part of tan$, and £ £ (0,1) can be written as a 
continued fraction 



£ = l/(ai + l/o2 + •••)) = [ai,02,. 



(4.2) 



For irrational £ the sequence [01,02, ■ ■ ■] is infinite and its finite approxi- 
mate is 

£„ = [01, . . . , a„] = p„/q„ (4.3) 

with mutually simple p n , q n . In the following, we use three important results 
on continued fractions ( Levy [1937| ) : 



(a) Estimate for the accuracy of the n-th convergent 

\€-Pn/Qn\ < l/?n?n+l , 



(4.4) 



(b) Asymptotic property of the sequence {aj} 



lim(o 1 ...a„) 1/ " = n( 1 



fe=i 



1 



k 2 + 2k 



lnfc/ln2 



= 2.63. 



(4.5) 



(c) Asymptotic property of the sequence {qj} 



lim -lng„ = 7T712 In 2 = 1.186. 



(4.6) 



Formulas (4.5) and (4.6) indicate a remarkable scaling structure of con- 
tinued fractions 

01 . . .o„ ~ A"g a (n) 



(4.7) 
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with scaling parameters A a , X q and some subexponential functions g a {n), g q (n). 

For an arbitrary irrational trajectory with an angle consider its n-th 
approximate and a corresponding rational trajectory with an angle 
such that tam9(") — ao + £„. This trajectory is periodic with the period 
Galperin and Zcmlyakov [199(| 



T n = const • q n (4.8) 

with the const = 1 or 2. The approximate trajectory is close to the original 
one for a time much larger than T n due to the fast convergence (4.4). That 
is just the property that one can see from the samples in FIGURE 2.5. The 



property (4.8) was used in Zaslavsky and Edclman [2001] to get the scaling 



for different quasi-periods of any irrational trajectory 

T n ~ X^grin) (4.9) 

which is similar to (4.7) with the time-scaling parameter At 

In A T = In A, = tt 2 /12 In 2 = 1.186 . . . (4.10) 

Simulation for i3 = 4.153087. . . gives In At = 1.17 ± 0.08 in an excellent 
agreement with (4.6). 

Expression (4.10) defines hierarchical structures of different flights in 
time. Corresponding hierarchical structures should exist for the lengths £ n 
of the flights along y. Nevertheless, there exist an ambiguity of the flight 
lengths l n since trajectories can have strong coordinate oscillations for the 
same time duration of flights. To find a corresponding scaling parameter 
X y , we assume that 



In A y = In Aden (4.11) 

where Aden means scaling parameter of different possible denominators of 
the rational convergent of £, obtained at the same hierarchical level along 
a trajectory, and the bar means averaging over such possibilities. For ex- 
ample, at the hierarchical level n, t; n has q n as the minimal denominator 
and ai . . . a n as the maximal denominator since 

oi . . . o„ > q n ■ (4-12) 

Applyiung (4.7) we obtain 

minA d cn = AT, maxA dcn = A Q (4-13) 
The simplest estimate is 

1 



In A, = In Aden « -(In A r + In A„) » 1.07... (4.14) 

The obtained information can be applied to construct the kinetic evolution 
of an ensemble of trajectories. 

All the numbers (4. 10), (4. 14) should be the same for the billiards with 
a slit (FIGURES 2.2(a), 2.3(a)) and with equal squares (FIGURES 2.2(b), 
2.3(c)). 
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FRACTIONAL KINETICS IN 
BILLIARDS 



The fractional space-time kinetic equation was introduced in Zaslavsky 
1992| and detailed fo r dynamical systems in Zaslavsky 1994a , EJ7 Saichev 



and Zaslavsky [1997] in order to describe self-similar non-Gaussian pro- 



cesses with strong intcrmittcncy. For a billiard of the type in FIGURES 
2.2(a,b) and 2.3(a-c), the equation has the form 



dtP 



V- 



d\y\ a 



(5.1) 



where f(y,t) is a distribution function, and fractional derivatives can be 
considered using their Fourier transform (F.T.): 



F.T. 



F.T. 



d\yf 



Usually, we are interested in asymptotics t — > 
moments 

(M 2m } = / y 2m f(y,t)dy 



= i\k\ (5.2) 
oo, \y\ — > oo. Since the 
(5.3) 



diverge for 2m > a ( Zaslavsky [19944 b| , |Saichcv and Zaslavsky [1997| ) , 
it is convenient to consider truncated distribuiton f(y,t), (y < y max ) and 
moments with m < maxm and maxm > 2, all of which are finite and 
present a self-similar evolution (see for example FIGURE 3.2 which shows 
m = 1, 2, 3, 4). With this comment we can get from (5.1) 



and expect that 



with 



(|y| Q ) = const, t 

(\y\ 2m ) ~ *" (m) 
fi(m) ~ m/i(l) 



P 



(5.4) 

(5.5) 
(5.6) 



(compare to (3.2)). 

In fact, the situation is more complicated ( ^aslavsky and Edclman [2000 |). 
The self-similarity of f(y,t) means that Eq. (5.1) is invariant under the 
transformation 

t-*\ T t, y -> \ y y (5.7) 

with appropriate values of the scaling parameters At, A. y . Then from (5.1), (5. 7) 
we obtain the fix-point condition 



lim {\ a y /4r 



1 



(5. 
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or 

(3 In A T = a In \ y + 2nik (5.9) 

with an integer k. Then the solution for (5.1) can be written in the form 
flZaslavsky and Edclman [2001[ |, Zaslavsky and Edclman [2U()(| ) 

(\y\ a )= £ C fc i* (5.10) 

k— — oo 

p k = an(l)/2 + 2mk/\n\ T 

with 

/i = 2 In A y / In A T , (5.11) 

and some expansion coefficients The expression (5.10) can be trans- 
formed into 

oo / 1 / \ 

(\y\ a ) 2/a = #* E Vk cos 2 ^rr~ + ^ ( 5 - 12 > 



k=0 



with new coefficients 2\ and phases ipk- The last expression shows the 
so-called log-periodicity with a period 

r log = hiA T (5.13) 

and with corresponding terms T>k=jto that typically are small. For the self- 
similar behavior of the moments of f(y, t) we also expect 

{y2m) „ t ,(m) J- f) k cos 2nk^- + i, k ) (5.14) 

with coefficients T>u and phases tj)j-, and 

fi(m) ~ m/i(l) = m/i = 2 In Ay/ In At ■ (5.15) 

Expressions (5. 13), (5. 14) are just the ones that should be tested by sim- 
ulations. As it was mentioned in the previous section, the results should be 
the same for the billiards in FIGURES 2.2(a) (2.3(a)) and 2.2(b) (2.3(c)). 
The results of simulations are presented in FIGURE 5.1. In FIGURE 5.1(a) 
we observe power law decay of the distribution of Poincare recurrences with 
7 w 2.7. There are also some oscillations which make deviations from this 
value of 7. FIGURE 5.1(b) shows moments for m = 1,2,3,4 and the self- 
similarity with /i = /Lt(l) ~ 1.7. There are small deviations from the linear 
dependence of /u(m) on m for m = 2, 3 and 4, /j,(2) = 3.5; /i(3) = 5.5; 
M (4) = 7.4. 

Let us consider the formula (5.15) and use for At and X y the values 
predicted by the continued fraction theory (4.10) and (4.14). It gives \i = 
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FIGURE 5.1. Statistical properties of the square billiards: (a) density distribution 
function of the recurrences P(t); (b) moments of x (full lines) and y (point-dash 
lines) for m — 1,2, 3,4 starting from the bottom; (c) fluctuations of the second 
moment of x vs. time; and (d) their Fourier spectrum. The data are obtained for 
4,048 trajectories during the time 10 s for each. 



/Lt(l) = 1.8 which is in a good agreement with the results of simulation in 
FIGURE 5.1(b). The values of 

7»2.7, (isil.7 (5.16) 

are also in agreement with (2.10). It is more delicate issue to evaluate 
the log-periodicity. There are small oscillations in the second moments 
(x 2 ), (y 2 ) time dependence. In FIGURE 5.1(c) we show their amplitude 

A(x 2 ) = (x 2 ) - const.*'' . (5.17) 

These oscillations do not look absolutely random due to the presence of 
structures. To find it, consider the Fourier transform of (5.17), i.e. 

= f dt e 2 " K log »> 4 A{x 2 ) (5.18) 

FIGURE 5.1(d) shows the almost continuous spectrum with some regu- 
lar peaks. The characteristic width of the regular part of the spectrum is 
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Ak x ~ 2. From (5.13) and (4.10), the width of the spectrum in the decimal 
logarithm basis should be 1/ log 10 Xt ~ 2 in the excellent agreement with 
the data of simulations. Similar agreement was obtained in Zaslavsky and 



Edclman [2001] for the billiard of FIGURE 2.2(a) (2.3(a)) type. 



6 RHOMBIC BILLIARD 



This type of billiard is shown in FIGURE 2.2(c,d). Its one-dimensional pe- 
riodic continuation (FIGURE 2.4) can be interesting for different applica- 
tions such as wave/ray propagation in nonuniform media. Two-dimensional 
periodic rhombic scatterers with an external held and small values of the 



Lyapunov exponents were considered in Lcpri, Rondoni. and Bcncttin [2000 1 




log 10 t 



FIGURE 6.1. Moments and distribution density P(t) of the recurrences for the 
rhombic billiard. The data are obtained from 512 trajectories during the time 
10 10 for each. 



The most interesting case is the irrational billiard, i.e. the case with an ir 
rational ratio of the rhombic diagonals. The rhombic billiard should inherit 
properties of the triangular billiards considered in Artuso, Guarneri, and 



[Rebuzzini [2000 1, Artuso, Casati, and Guarneri [1997 [, Casati and Prosen 
|[1999 [. In FIGURE 6.1 we presented the recurrences distribution P(t) for 
which the slope recurrence exponent is 7 = 2 ± 0.1, and the coordinates 
moments (x 2m ), (y 2m ) versus time. The latter shows two exponents: dur- 
ing the time interval < t < 10 7 , /i(l) ~ 1.8 -j- 1.9, and after that time 
10 7 < t < 10 10 , fi(l) « 0.95±0.05. We cannot guarantee that this exponent 
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is not an intermediate one, but it satisfies the condition (2.10) and both 
values of 7 and fx(l) for large t are close to the values obtained in |Casati 



and Prosen [1999 1 for the triangular billiard. These values will be discussed 



in Section 8. 




Ixl » 10" X 



FIGURE 6.2. Rhombic billiard with the x-diagonal -=1.6 and {/-diagonal = 
1.2290472 (a) a sample trajectory in a square shows slow evolution; (b) phase 
plane £,ve (I is the upper side of the square in (a)) shows slow evolution in the 
angle (i.e. along vi); (c) a sample trajectory in large scale; and (d) its zoom. 



In FIGURE 6.2 we present samples of trajectories which show trap- 
pings and very slow escape from the trapping domains. This feature of the 
trajectories can be especially well observed from FIGURE 6.2(b) of the 
Poincare section: each point corresponds to a trajectory crossing of the in- 
terval I : y = 1, x e (0, 1) in FIGURE 6.2(a). There is extremely slow 
process of mixing along the velocity vi and the filling of the phase space is 
performed mainly along x for some special values of vg. We will speculate 
on the recurrences estimate in Seciton 8, applying for this type of dynamics. 
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7 BACK TO THE CONTINUED 
FRACTIONS 

Here we present a comment to the continued fractions properties. Since 
trajectories in the billiard with a slit or a square bear scaling properties of 
the continued fractions from one side, and the log-periodicity from another 
side, one can expect that the statistical features of continued fractions 
should have similar type of the log-periodic oscillations. 

Consider a large number (ensemble) of irrational numbers < 1}. For 
each representative £^ of the ensemble, consider their n-th approximate, 

^ = [a^,a^,...,a^]=p^/ q ^ . (7.1) 

The new ensemble is the ensemble of the denominators {q^}. We can 
introduce a distribution function &(q; n), i.e. a probability density to have 
value q for the denominator of the convergent of irrational numbers at n-th 
step. The variable n plays a role of the discrete time, <I>(q; n) should be a 
universal function for which we can introduce moments 

(ll m ) = J q 2m mn)dq (7.2) 

The asymptotic property (4.7) suggests that 

±-\n(q^)=nc(l + d n ) (7.3) 

with 

c = lnA 9 (7.4) 

and d n depends slowly on n. The conjecture is that the Fourier spectrum 
of d n consists of two parts: quasi-periodic and continuous, similar to what 
we have for the billiards. 

To verify the conjecture, we considered 10 3 different irrational numbers 
£0) < i (y — 1, . . . , 10 3 ) and their approximates with the quadrupole 
precision. The behavior of moments (q 2m ) is presented in FIGURE 7.1 
up to n = 150. The plot shows oscillations which are stronger for higher 
moments. The slope corresponds to the value (7.4) with an accuracy up to 
the level of oscillations. The value of oscillations with respect to the straight 
line is shown in FIGURE 7.1(a), and their Fourier transform is in FIGURE 
7.1(b). The figures show a finite, quickly decaying peaks of modulation of 



the continuous type spectrum. The last one has been proved in tbragimov 
and Linnik [197l| ]. 
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FIGURE 7.1. Statistical features of the continued fractions: moments of the de- 
nominators {q„ m ) vs. number n of the approximation level and fluctuations of 
the 8-th moment with the spectral function. 



8 ESTIMATES FOR THE 

DISTRIBUTION OF RECURRENCES 

In this section we provide rough estimates for the distribution of recurrences 
P(t) or, more accurately, for possible values of the recurrence exponent 7 
when the P(t) behaves algebraically. As it was mentioned in (2.6) the value 
of 7 should exceed 2 for the bounded Hamiltonian dynamics. 

A qualitative way to obtain P(t) is to consider an enveloped phase volume 
r(i) evolution with time and a probability of a particle to return back from 
this volume to the initial one To < T(t). One can consider a small ball 
of the size e that moves in phase space and cover it with time. There are 
different possibilities: 

Fi(i) = const. • t 1+Sl 
f 2 (t) = const. • t 3/2±S2 
f 3 (t) = const. • t 2±Sa 

f 4 (t) = const. • t 3±5i (8.1) 

where Sj are considered as corrections to four different main values. The 
case of T\(t) corresponds to almost one-dimensional coarse-grained dynam- 
ics. It seems that this case happens in the rhombic and triangular billiards 
and FIGURES 6.2(b,c) confirm this statement. The same case can appear 



also in one-dimensional dynamics (sec Young [199S ] and some rigorous re- 
sults therein). 

The case of f^f) corresponds to the linear evolution along one coordi- 
nate and the diffusional type evolution along the other coordinate along 
which the phase volume growth is proportional to i 1 / 2 . This is one of the 
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most typical case for many systems with intermittent chaotic dynamics 
flZaslavsky and Edclman [2000] ). The third case of T^^t) corresponds to 
the almost linear phase space coverage along one coordinate and then the 
similar coverage along another coordinate. The alternations of the coordi- 
nates can be randomly distributed but the result provides almost £ 2 -growth 
of Ts(t). It seems that it is the case of the original Sinai billiard when 63=0 
( phernov and Young [2000 1). The case of f 4(f) can appear in speci al situa 



tions of 1 1/2 degrees of freedom or in higherc dimension system (Rakhlin 



[2001| ). 

For (8.1) we obtain the return probability during the time interval (0,t) 
as To/Tj(i) and finally 

P^A^- const./*- (8.2) 

with the corresponding values of jj 



(8.3) 



where Sj represents the corrections to the main exponents. We have to 
recall that for some cases exact values of the exponents have no meaning 
since the log-periodic oscillations of the distributions P(t). 




9 CONCLUSIONS 

The randomness with zero Lyapunov exponent, that we call pseudochaos, 
can appear in numerous applications. Its analysis provides a link between 
the structure of streamlines of flows, magnetic field lines, billiards, con- 
tinued fractions, and fractional kinetics. It seems that the dynamical pro- 
cesses that can be described by the fractional type kinetic equation reveal 
a polynomial complexity, while the chaotic hyperbolic dynamics has the 
exponential growth of the coarse-grained phase volume. 
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